Evolution Letters

Does the distribution of fitness effects of new mutations vary across the genome? Under the classical Fisher Geometric Model (FGM) we might not expect it to. In FGM, phenotypic traits are envisioned as dimensions of a landscape, with fitness determined by position in the landscape, i.e., the particular combination of traits of an individual. New mutations are represented by vectors that move from an ancestral to a new phenotype. In classical FGM these vectors affect all trait dimensions simultaneously (universal pleiotropy). However, introducing partial and modular pleiotropy into an FGM framework leads to an expectation that parameters of the DFE will vary with mutational pleiotropy-the number of traits affected by individual mutations. Here we address this prediction by investigating whether traits related to mutational pleiotropy, expression level and network connectivity, affect the parameters of the DFE using whole genome data from A. thaliana and C. grandiflora, two closely related Brassica species that vary significantly in their demography and mating system, and therefore, in effective population size and the effects of linked selection. Results were similar across both species. We found that expression level and network connectivity were predictive of the parameters of the deleterious DFE, even once co-correlations among genome biology traits were accounted for. Our results suggest that, across the genome, molecular evolutio(high mutational pleiotropy). nary patterns agree with the predictions of FGM, albeit relaxing the assumption of universal pleiotropy, and that variation in mutational pleiotropy among genes is sufficient to have detectible effects on the DFE. Significance statementHow do the effects of new mutations vary across the genome? If mutations in some genes affect many traits (high mutational pleiotropy), we hypothesise they will be more strongly deleterious, with lower variance in their selective effects. We test this by investigating the distribution of effects of new mutations across genes that vary in features that are related to mutational pleiotropy: expression level, gene network connectivity, and number of associated GO terms. The mean strength and coefficient of variation of selection of new mutations varied across genes with different features in the manner expected by our hypothesis. This demonstrates that important parameters of molecular evolution can vary across the genome with genome architecture.

While life-history strategies are typically fixed within species, evolutionary transitions between perenniality and annuality can occur. In clonal seagrasses, annual and perennial plants often coexist in the same population, providing a unique model for studying the genetic basis of this transition. Two seagrass Zostera marina populations in South Korea display a striking dichotomy: shallow-water sub-populations follow a typical perennial strategy, whereas their deep-water counterparts are annual. Here we show that this shift from perenniality to annuality, potentially caused by the SAPK7 gene, is genetically coupled with the CAO gene, which is under strong positive selection for low-light adaptation. The up-regulation of the SAPK7 gene triggers early flowering in seedlings, before the formation of any lateral shoots via asexual reproduction. In this special case where the genet contains only one ramet, the post-reproductive death of the ramet is equivalent to the death of the whole genet, which explains the annual phenotype. Our findings reveal a mechanistic example where annuality overcomes perenniality by hitchhiking on a positively selected gene. Given that the ancestral state of plants is perennial, this coupling of annuality with beneficial alleles may represent one of the pathways for the repeated evolution of annual life histories across flowering plants.

The potential for environmental change to compound the detrimental effects of inbreeding depression in small and isolated populations is a significant concern in conservation biology. Previous evidence syntheses suggested that environmental stress exacerbates inbreeding depression, but were based on limited data. Here, we comprehensively test the relationship between inbreeding depression and environmental stress in natural populations using Bayesian mixed-effects meta-analysis on a large, high-quality data set of 2127 inbreeding depression effect sizes from animals and plants. Our results show that inbreeding depression is significantly higher in benign than in stressful environments. Analyses of both inbreeding depression and stress-induced changes in genetic load supported a unimodal (humped) relationship between the costs of inbreeding and stress intensity, with a peak at intermediate stress. At the highest levels of stress there was, on average, a significantly greater inbreeding load in benign than in stressful environments. We suggest that the lower cost of inbreeding associated with extreme stress results from constraints on the expression of inbreeding depression as fitness and phenotypes decline towards zero. Our findings help to resolve long-standing uncertainty around how inbreeding and environmental change interact, revealing that inbreeding responses vary non-linearly with environmental stress intensity, but showing that stress does not generally amplify inbreeding depression. As such, they will inform both the management of populations of conservation concern and predictions of species responses to global environmental change.

Antibiotics impose strong selection on bacteria, resulting in the emergence and spread of antimicrobial resistance. To date, the consequences of resistance evolution for other traits, especially virulence as a key life-history characteristic of high medical relevance remain poorly understood. A central limitation of existing work is reliance on clinical isolates with complex evolutionary histories, hindering causal inferences on resistance-virulence relationships. Such causal information is critical for our understanding of the distribution of resistance-associated evolutionary trade-offs, which may additionally guide optimization of treatment designs. The objectives of our study are to address these knowledge gaps using an experimental evolution framework with the human pathogen Pseudomonas aeruginosa. We quantified changes in resistance, life-history characteristics, in-vivo virulence with a Caenorhabditis elegans infection model, and whole genome sequences for bacterial clones, which had been independently evolved under three antibiotics with distinct cellular targets. We found that the resistance-virulence trade-off depended on the used antibiotic and was additionally driven by the underlying evolutionary path to resistance. Evolved resistance to the fluoroquinolone ciprofloxacin correlated positively with virulence, dependent on genetic changes in either the antibiotic target or efflux regulation. Conversely, evolved resistance to piperacillin/tazobactam, a {beta}-lactam/{beta}-lactamase inhibitor combination, was negatively correlated with virulence, contingent on the coincidental spread of resistance mutations and a large genomic deletion, containing numerous virulence genes. Lastly, evolved resistance to the aminoglycoside streptomycin led to only minor virulence changes. Overall, these antibiotic-specific evolutionary trajectories challenge the assumption of a universal resistance-virulence trade-off and demonstrate that antibiotic choice itself shapes pathogen virulence potential. Significance StatementAntimicrobial resistance and bacterial virulence are typically studied in isolation, yet they are shaped by the same evolutionary pressures and may directly compete for cellular resources. We addressed a fundamental but unresolved question: does evolution of resistance to an antibiotic make a pathogen more or less dangerous to a host? By characterizing experimentally evolved Pseudomonas aeruginosa, a leading cause of drug-resistant infections, we show that the answer depends on which antibiotic drove resistance. A positive resistance-virulence relationship occurred upon adaptation to a fluoroquinolone antibiotic, whereas it was negative upon {beta}-lactam resistance, and unchanged upon aminoglycoside resistance evolution. These antibiotic-specific evolutionary trajectories challenge the assumption of a universal resistance-virulence trade-off and suggest that antibiotic choice shapes pathogen danger beyond drug susceptibility.

Gene regulatory divergence has emerged as a key feature in speciation, influencing gene expression differences that accumulate between diverging populations. Transcriptional regulation, mediated by cis- and trans-acting factors, modulates diverse developmental processes and is responsible for distinct species-specific gene expression profiles. Within interspecies hybrid individuals, negative interactions between divergent cis- and trans-acting factors can lead to gene misregulation and hybrid dysfunction at the organismal level. Such gene regulatory mismatch might disproportionately impact sex-biased and tissue-biased gene regulatory networks due to their unique selective pressures. To address these issues, we investigated the role of regulatory divergence in asymmetric hybrid incompatibility between sister species of Caenorhabditis nematodes (C. remanei, C. latens) by analyzing gene expression of reciprocal hybrids for each sex and key tissue types. Despite severe hybrid male sterility, hybrid males showed less misexpression of sex-biased genes than hybrid females, suggesting that the organismal phenotypic outputs of male-biased gene regulatory networks are more vulnerable to disruption than female-biased genetic networks. Additionally, we found more genes associated with cis- than trans-regulatory divergence, supporting the notion of a disproportionate role for cis-regulatory divergence between species. Moreover, we document extensive cis-trans compensatory X-linked regulatory divergence specifically from male transcriptomes, indicating distinct molecular evolutionary outcomes of stabilizing selection on regulatory controls in males and females. These insights derived from asymmetric hybrid misexpression expand our understanding of the evolution of sex-biased gene regulation in the face of stabilizing selection and identify candidate genes contributing to Caenorhabditis post-zygotic reproductive isolation.

Bird migration evolved as an adaptation to seasonally changing habitats. Migratory behaviour can vary within the same species in case of partial migratory behaviour, i.e. one population (or individual) is migratory and another one is resident. Species that exhibit a wide variety of migratory phenotypes provide valuable systems to understand the evolutionary drivers behind different phenotypes and how populations adapt to habitats with distinct seasonality. The European robin (Erithacus rubecula) expresses migratory behaviour in central and northern areas of the species distribution range, whereas populations in the South and on the Macaronesian islands are predominantly resident, providing a suitable system to investigate these questions. We use high coverage whole genome re-sequencing data of 125 European robins to investigate how migration behaviour affects population structure and demography, and how it affects the selection landscape in the genome. Genetic structure in European robins coincides with migratory phenotype and geography and populations are characterised by distinct demographic histories. Our results suggest that both the continental resident population as well as the Macaronesian island populations have derived independently from an ancestral migratory population. Unexpectedly, tests for differential selection revealed extensive positive selection pressure acting across all chromosomes in the resident populations, while selective sweeps are largely absent from migrants. We speculate that this might be an analytical artifact due to mismatching timescales between what population genomics methods can detect and the scale on which migration behaviour likely evolved in the robin. We suggest that future studies on the genomics of migration should more focally account for different time scales on which these processes happen, such as including the wider phylogenomic background of the target species, to capture the full evolutionary history of migratory traits.

Empirical studies have demonstrated that mutation rates may change with individual condition, such as in the case of stress-induced mutagenesis. This has led to the hypothesis that condition-dependent (or "plastic") mutation rates could be selectively favoured, as the increased production of new mutants in times of maladaptation enhances evolvability, the ability to undergo adaptive evolution. However, while empirical evidence for condition-dependent mutation rates is accumulating, theoretical models studying their evolution are lacking. Here, we employ an individual-based simulation approach to examine the evolution of condition-dependent mutation rates in a changing environment. We find that condition-dependent mutation rates consistently evolve when the environment changes at an intermediate pace. Furthermore, populations with condition-dependent mutation rates are substantially better adapted to their (changing) environment. Finally, the evolutionary dynamics of condition-dependent mutation rates are both accelerated and destabilised when the mutation rate is self-referential (i.e., when mutator loci affect their own mutation rate). We conclude that condition-dependent mutation rates (and thus evolvability) can readily evolve in changing environments. Significance statementMutation provides the raw material for evolution. Mutation rates thus tune evolvability, the ability to undergo adaptive evolution: if mutation rates are too low, evolution is impeded; if mutation rates are too high, adaptive traits cannot be maintained. Using a theoretical model, we explore the evolution of plastic mutation rates that systematically depend on the condition of the organism and its environment. An example is stress-induced mutagenesis in bacteria, which is implicated in the evolution of antibiotic resistance. We show that plastic mutation rates readily evolve, providing "well-timed" variation specifically when organisms are poorly adapted. Such plastic mutation rates thus facilitate better adaptation to changing environments, and their evolution provides an example of evolvability itself evolving.

If an individuals niche is determined by its genotype, then competition for limiting resources should be most intense among individuals of the same genotype. Theory predicts this will act to maintain genetic variation, but whether this mechanism operates under natural conditions is unclear. Using long-term observations of a population of free-living Trinidadian guppies, we asked (i) whether competition was strongest between kin, and (ii) whether this process maintained genetic variation. Competition between kin was 1.5-1.8 times stronger than that between non-kin. This contributed to balancing selection: after [~]10 generations, variation was 29% higher than expected under drift. Our results show that relatedness can play a major role in structuring ecological competition, with broader consequences for theories of inclusive fitness. One-sentence summaryHeritable variation is maintained due to resource competition being more intense among kin.

Climate change can threaten population viability by disrupting sex ratios in species whose sex is influenced by temperature. While species with sex chromosomes were historically considered immune, in some species, temperatures can override genetic sex determination via sex reversal, leaving them vulnerable to climate-driven sex ratio shifts. The Tasmanian spotted snow skink (Carinacincus ocellatus), a viviparous reptile with an XX/XY system, provides a compelling case study. While laboratory studies demonstrated that extreme thermal conditions induce female-to-male sex reversal (XX males), its occurrence in the wild remains unexplored, limiting our understanding of actual climate impacts. Integrating 23 years of phenotypic and genetic sexing data across two climatically distinct populations, we provide the first evidence of sex reversal in a wild viviparous reptile. XX reversal occurred in both populations, affecting up to 23.5% of XX births in the warmer population, and was associated with colder minimum daily temperatures. Despite high birth rates in some years, sex-reversed adults were rare. We also identified putative XY females, suggesting bidirectional sex reversal and reinforcing the extreme plasticity of reptilian sex determination. Ultimately, sex reversal could act as an evolutionary trap, potentially compromising population viability as climate instability increases.

Hosts defend themselves against parasites through resistance (reducing parasite burden) and tolerance (reducing the fitness cost of infection without affecting parasites). This distinction has important evolutionary implications: resistance is predicted to maintain polymorphism while tolerance tends to fix, and only resistance is expected to provoke parasite counter-adaptation. The reaction-norm framework, which infers tolerance from the slope of host fitness regressed on parasite burden, assumes that a shallow slope reflects parasite-independent host protection. We test this assumption using a within-host model in two variants: microparasites (Model 1, with within-host replication) and macroparasites (Model 2, without). Sublethal immunity impairs the host-exploitation rate of the parasite, reducing both growth and per-parasite virulence without killing them. We show that this generates systematic slope differences among host genotypes that the framework interprets as variation in tolerance. Furthermore, the ranking of slopes across genotypes reverses between linear and sigmoidal damage functions: under linear damage, the strongest immune responder appears most tolerant; under sigmoidal damage, the weakest responder does. Decomposition of the damage reduction shows that virulence reduction accounts for the majority of the effect across both model variants. Thus, the reaction-norm slope cannot determine whether host fitness is maintained by parasite-independent tissue protection or by sublethal impairment of parasites.

Many proposed mechanisms for the evolution of cooperation among unrelated individuals rely on relatively demanding cognitive abilities that are not widespread across taxa. In contrast, individual heterogeneity is a pervasive feature of animal groups, encompassing differences in personality as well as physical and cognitive traits. Such heterogeneity can promote the evolution of cooperation, yet its role has received comparatively little attention, particularly as a source of variation giving rise to social organization such as leadership. A specific form of leadership can emerge under unstable environmental conditions, when some individuals become better suited than others to initiate action and influence the behavior of their peers. Unlike fixed dominance hierarchies, emergent leadership can rapidly adjust to changing environmental conditions, thereby reshaping group organization. Because it does not require the maintenance of stable hierarchies, this form of leadership can arise even in species that do not have the cognitive capabilities to sustain complex social structures. In this work, we investigate the combined effects of individual heterogeneity and emergent leadership on the evolution of cooperation using an evolutionary game-theoretic model in which individuals may assume the roles of leaders or followers according to their strength, representing individual differences in suitability to prevailing environmental conditions. We examine different levels of population heterogeneity together with increasingly complex strategy sets requiring progressively greater informational requirements, allowing individuals to condition cooperation on their own strength, leadership role, or both. Our results show that the interplay between leadership and heterogeneity promotes the evolution of cooperation, particularly when only a small fraction of individuals act as leaders. Under these circumstances, cooperation evolves even when individuals employ the simplest possible strategies. Under harsher ecological conditions, cooperation can be sustained by more sophisticated strategies, specifically by conditional strategies that prescribe cooperation when individuals are strong or leading and defect when acting independently. Author summaryIn this study, we propose that emergent leadership mediated by individual diversity can boost the evolution of cooperation in animal groups. Building on growing evidence on the heterogeneity of animal capabilities and personalities, we focus on the fleeting leadership that emerges in animal groups when facing rapidly changing environmental conditions. We suggest that this type of leadership that emerges from individual differences in strength--a generic quality encompassing those characteristics that make an individual more fit to lead in a given situation--does not require complex cognitive capabilities from the animals and represents a valid alternative to more demanding strategies proposed in the past to explain the evolution of cooperation. Using an evolutionary game theory model, we show that if a population includes a few strong players, these can become influential leaders and guide the actions of their peers to achieve cooperation. Although the naive strategy of always cooperating is sufficient for cooperation to evolve, the introduction of more complex strategies leads players to cooperate only when they are more likely to be recognized as influential leaders. These strategies are more effective in promoting cooperation under unfavorable ecological conditions and are also more robust against exploitation by defectors.

Conjugative plasmids drive the dissemination of antimicrobial resistance genes and other conditionally beneficial accessory genes, but otherwise inflict fitness costs on their hosts that limit their spread. These costs can be ameliorated by compensatory mutations on the chromosome, the plasmid, or both combined (combined compensation). Mutants with plasmid-borne or chromosomal compensation differ in how they spread and compete, making the location of compensation an important determinant of plasmid and antimicrobial resistance (AMR) dynamics. We collate experimental data on plasmid-host co-evolution which, albeit limited, suggest compensatory benefits differ by location and are highest for combined compensation. We develop and analyse a mathematical model of compensatory evolution, finding that the long-term location of compensation is mainly determined by the highest cost reduction. Short term, however, succession dynamics arise from differences between locations: for instance, plasmid-borne compensation spreads horizontally, initially dominates, and can even facilitate the establishment of chromosomal or combined compensation. Strong trade-offs between compensation and either resistance or conjugation render compensation non-viable, but only conjugation trade-offs are location-dependent, disadvantaging plasmid-borne compensation and generating oscillatory dynamics. Our findings suggest plasmid-borne compensation may initially accelerate AMR-plasmid spread, whereas long-term chromosomal or combined compensation may enable hosts to accumulate multiple AMR-plasmids, promoting multidrug resistance.

Introgression following hybridization is increasingly recognized as a major driver of evolution. However, its importance depends on its frequency in nature, which remains to be quantified. To address this, we provide a snapshot of ongoing introgression in a whole species assemblage (4126 ant colonies). 23% of all 82 local species show signs of introgression, which is more than twice previous estimates. Introgression is typically subtle, yet we find that it contributes measurably to genetic diversity. Species divergence, rather than classical prezygotic reproductive barriers (mating phenology, ecological niche, fine-scale habitat use) constrains introgression, suggesting that the main reproductive barriers are postzygotic at this stage of divergence. Our results indicate that introgression may be a common but often overlooked feature of natural communities.

While asexual species can often outcompete their sexual counterparts over ecological timescales, their long-term evolutionary success is hindered by a diminished ability to purge deleterious mutations and to adapt to changing environments. However, some asexual species persist for millions of years, and a major question in evolutionary biology is how they do so. One solution is to occasionally reproduce sexually, as has been shown in a handful of primarily asexual species. Here, we investigate the possibility of rare sex in the clonal raider ant, Ooceraea biroi. We report the whole-genome sequence of a previously uncharacterized clonal line and, using population genetic and phylogenetic analyses, show that it originated through sexual reproduction between two extensively studied clonal lines. The mitochondrial genome of this clonal line differs from that of the maternal clonal line at only a single nucleotide, suggesting that the sexual reproduction event occurred within the past few hundred years. These results demonstrate that sex occurs sporadically in the clonal raider ant, allowing it to generate new genetic combinations and potentially to overcome some of the costs of asexuality.

Biological invasions on islands provide a natural framework to study how dispersal and connectivity influence evolutionary and ecological processes. The avian nest parasitic fly, Philornis downsi - first recorded in Darwin's finch nests in 1997 - causes high mortality in endemic land birds, yet its inter-island and sex-specific patterns of dispersal and genetic structure remain poorly understood. We use low-coverage whole genome sequencing to investigate genome-wide patterns of genetic diversity, directional migration and effective population size in P. downsi across five major Galapagos Islands and its native range in mainland Ecuador. We find evidence for a genetic bottleneck in the Galapagos, isolation by distance, and evidence that the island closest to the Ecuadorian mainland, San Cristobal, is genetically divergent from the other four islands sampled, despite retaining the highest genetic diversity. No evidence was found for sex-biased dispersal; however, sex-biased genetic structure was detected using only markers from inferred autosomal scaffolds. We found asymmetric gene flow with higher migration rates from San Cristobal westward to the other islands, matching the direction of both southeast trade winds and major cargo shipping routes. Our results suggest both natural and human-mediated colonisation of P. downsi from the mainland through San Cristobal to the other islands, followed by high inter-island dispersal among closely situated sink islands. Our findings are critical for prioritising islands for control strategies that will reduce P. downsi impacts on vulnerable endemic birds and underscore the value of understanding directional migration patterns for managing invasive species in metapopulations.

Sexual conflict has been hypothesized as a driver of speciation, though its effects are likely heterogeneous across phylogenies and between sexes. The semi-aquatic bug, which inhabits water surfaces across diverse aquatic environments, has long served as a model for studying sexual conflict. While previous studies have focused on rapid antagonistic coevolution and the genetic basis of sexually antagonistic traits, the macroevolutionary consequences of asymmetrical sexual conflict--particularly male-dominated grasping traits versus female resistance--remain largely unexplored. Within the subgenus Pseudovelia, males exhibit pronounced phenotypic diversification in grasping structures, whereas females show modest, clade-specific resistance traits, suggesting male-biased asymmetric conflict. This system presents a valuable opportunity to examine how sexual conflict influences diversification and asymmetrical trait evolution across lineages. Using 204 individuals, representing over half of the subgenus's species diversity, we reconstructed a time-calibrated phylogeny, quantified diversification rates, assessed sexual conflict intensity across clades, and analyzed correlations between sexual trait evolution and diversification. Our results reveal extensive phylogenetic conflict, particularly within the East Asian clade, driven by introgression and incomplete lineage sorting (ILS). Furthermore, we observe significant phylogenetic heterogeneity in both phenotypic evolution and diversification rates. Notably, a male "trait package" enhancing grasping ability likely drives rapid diversification in the recently radiated "South China" lineage. In contrast, grasping traits involving abdominal segment VIII are associated with lower conflict intensity, facilitating greater evolutionary flexibility in female resistance and resulting in lineage-specific counter-adaptations. These findings highlight the heterogeneous dynamics of asymmetrical sexual conflict in shaping diversification and speciation.

High investment in current reproduction can limit future reproductive opportunities, but how selection shapes these hormone-mediated traits remains poorly understood. Androgens can mediate male reproductive investment, and in three-spined stickleback (Gasterosteus aculeatus), exert antagonistic effects on breeding effort versus spermatogenesis. To understand how shifts in reproductive strategy shape this tradeoff, we compared testes transcriptomes and androgen production between two recently diverged stickleback ecotypes that differ in reproductive strategy: the ancestral "common" ecotype, which provides paternal care, and the non-parental "white" ecotype, which has lost paternal care and prioritizes mating effort. During typical breeding, testes gene expression differed little between ecotypes. However, under prolonged summer-like conditions, testes gene expression diverged substantially. Common-biased genes were enriched for meiotic functions and spermatogenic cell type markers, suggesting commons had initiated spermatogenesis while whites had not. Instead, whites expressed higher levels of steroidogenic candidate genes and released significantly more 11-ketotestosterone than commons, indicating sustained investment in current reproduction. F1 hybrids released 11-ketosterone at intermediate rates, suggesting a genetic basis for this divergence. Sustained androgen production in whites may possibly delay the transition into spermatogenesis, limiting investment in future reproduction. These results illustrate how selection on hormonally-integrated traits can drive rapid divergence in life history strategy.

Shifts in flowering phenology are one of the most well studied plant responses to global climate change. Many studies have documented these shifts and their drivers, including some that describe altered patterns of co-flowering among taxonomically broad species within communities. In comparison, few analyses have examined systematic changes in co-flowering between closely related, interfertile species, where co-flowering can have unique evolutionary consequences. To address such shifts in co-flowering among close relatives, we investigate phenological responses to climate change and its effect on patterns of co-flowering over the past 124 years in 52 species of North America violets (Viola). This genus has many co-occurring species that reproductively interact via shared pollinators and hybridization. We use ~14,000 herbarium records along with environmental and species trait data to model the magnitude of recent flowering phenology shifts, environmental variables and/or species traits associated with these shifts, and resulting changes in co-flowering among species. While both the magnitude and direction of phenological shifts varied among Viola species, nearly half (25/52) show significant changes in flowering day. Regardless of whether flowering was advanced or delayed, flowering date was most consistently associated with local mean temperature. Of six species-level traits, geographical region also significantly predicted flowering shifts, consistent with environment and geography together explaining broad phenological responses across this group. These shifts have produced significant changes to pairwise patterns of co-flowering among species -- ranging from a 59 day increase in co-flowering to complete loss of co-flowering overlap. Sympatric pairs specifically have experienced both increases and decreases in co-flowering, with a geographic pattern of increased co-flowering occurring mainly in eastern US and decreased co-flowering common in western US. Because Viola species are generalist pollinated and already known to hybridize, these new co-flowering patterns could further undermine reproductive barriers among species in this genus.

BackgroundLyme borreliosis is a common tick-borne disease in Europe caused by spirochetes of the Borrelia burgdorferi sensu lato complex, including Borrelia afzelii, which is maintained in nature through interactions with rodent reservoir hosts. These spirochetes have evolved several surface proteins to manipulate rodent host immunity, some of which remain polymorphic in Borrelia populations. Among these proteins, OspE, which binds the host complement-regulating factor CFH to evade destruction by complement, is one of the most variable. Yet, what evolutionary forces maintain this polymorphism is not well understood. Motivated by a recent discovery of CFH polymorphism in the bank vole (Clethrionomys glareolus), the main reservoir host of B. afzelii, we hypothesized that the polymorphism is maintained by host-parasite coevolution involving specific associations between host and parasite genetic variants. MethodsWe analyzed associations between bank vole CFH alleles and B. afzelii OspE variants across three datasets sampled in Poland. Selection acting on OspE was evaluated using omegaMap. Host-pathogen genotype associations were tested using partial redundancy analysis (RDA), and co-structure was assessed using co-correspondence analysis (CoCA). ResultsWe found that OspE evolves under positive selection, however, we found no evidence for an association between OspE and host CFH variants at the individual level based on RDA or at the population level based on CoCA. ConclusionsDespite evidence of positive selection acting on OspE, we found no support for specific genetic matching between B. afzelii and its bank vole host at the CFH-OspE interface. These results suggest that the evolution of CFH and OspE may be shaped by broader selective pressures, potentially including interactions with multiple host species.

Ancient whole-genome duplications (WGDs) are thought to have played a major role in plant evolution, but robust inference of the patterns and drivers of polyploid establishment through deep time remains challenging. We re-evaluate the recent large-scale study by Chen et al. (2026), which linked polyploid establishment throughout angiosperm evolution to periods of climatic instability and low species richness. We identify five conceptual and methodological issues that substantially affect these conclusions and collectively undermine the proposed temporal and ecological associations. We hope that clarifying these issues will support future efforts to understand the evolutionary role of polyploidy in plant evolution.